Half-metallic ferromagnetism induced by dynamic electron correlations in VAs 
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The electronic structure of the VAs compound in the zinc-blende structure is investigated using a 
combined density-functional and dynamical mean-field theory approach. Contrary to predictions of a 
ferromagnetic semiconducting ground state obtained by density-functional calculations, dynamical 
correlations induce a closing of the gap and produce a half-metallic ferromagnetic state. These 
results emphasize the importance of dynamic correlations in materials suitable for spintronics. 
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There is a growing interest in finding potential ma- 
terials with high spin polarization. Suitable candidates 
are Half-metallic ferromagnets (HMF) which are metals 
for one spin direction and semiconductors for the other 
[HQ As a result, HMF are expected to show a 100% 
polarization and can in principle conduct a fully spin- 
polarized current at low temperature. Therefore, HMF 
are ideal candidates for potential applications in spin- 
tronics 0, El • These materials are equally interesting for 
basic research, in particular concerning the origin of the 
half-metallic gap, the nature of inter-atomic exchange in- 
teractions Q, its stability at elevated temperatures H, 
or the effect of electron correlations Q. 

Most theoretical efforts for understanding HMF are 
supported by first-principles calculations, based on 
density- functional theory (DFT). In fact the very dis- 
covery of HMF was due to such calculations |l|. DFT 
calculations are usually based on the Local Spin Den- 
sity Approximation (LSDA) or the Generalized Gradient 
Approximation (GGA). These approximations have been 
proved very successful to interpret or even predict mate- 
rial properties in many cases, but they fail notably in the 
case of strongly-correlated electron systems. For such 
systems the so-called LSDA+f/ (or GGA+U) method is 
used to describe static correlations, whereas dynamical 
correlations can be approached within the LSDA+DMFT 
(Dynamical Mean-Field Theory) 0,0. An important dy- 
namical many-electron feature of half-metallic ferroma^ 
nets is the appearance of non-quasiparticle states 0,T 
which can contribute essentially to the tunneling trans- 
port in heterostructures containing HMF |10t lllj . 

Equally interesting materials for spintronics applica- 
tions are ferromagnetic semiconductors [l^ . . Candi- 
date systems are ordered compounds such as europium 
chalcogenides (e.g., EuO) and chromium spinels (e.g., 
CdCr2Se4) as well as diluted magnetic semiconduc- 
tors (e.g., Gai-^Mn^As) 0|. Unfortunately, all of them 
have Curie temperatures much lower than room temper- 
ature. On the other hand, VAs in the zinc-blende struc- 
ture is, according to density- functional calculations |l4j . 
a ferromagnetic semiconductor with a high Curie tem- 
perature. Unlike CrAs 0, CrSb 0, and MnAs [HI, 



VAs has not yet been experimentally fabricated in the 
zinc-blende structure, but the increasing experimental 
activity in the field of the (structurally metastable) zinc- 
blende ferromagnetic compounds is promising in this re- 
spect. 

In this Letter we investigate the effect of electronic in- 
teractions on VAs in the zinc-blende structure using dy- 
namical mean-field theory. Our main result is displayed 
in Fig. (especially the inset): While this material is ex- 
pected to be a ferromagnetic semiconductor from density- 
functional theory (LSDA/GGA) or static LSDA+J7 cal- 
culations, the inclusion of dynamic Coulomb correlations 
within the LSDA+DMFT approach predicts a metallic 
behavior, due to the closure of the gap in the majority- 
spin band. Moreover, since the minority-spin band gap 
remains finite, the material is found to be a half-metallic 
ferromagnet. To our knowledge, this is a first example in 
which dynamic correlations transform a semiconductor 
into a half metal. This remarkable result demonstrates 
the relevance of many-body effects for spintronic materi- 
als. 
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FIG. 1: Schematic representation of the p-d hybridization and 
bonding-antibonding splitting in VAs. 

The main features of the electronic structure of VAs 
[l4| are shown schematically in Fig. ^ The t2 g states 
hybridize with the neighboring As p states, forming wide 
bonding and antibonding hybrid bands. In contrast, the 
e g states form mainly non-bonding and narrow bands. 
The Fermi level falls between the e g and the antibonding 
t2 g in the majority-spin bands, and between the bond- 
ing t2g and the e g in the minority-spin bands. Thusc, 
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the material is a ferromagnetic semiconductor, showing 
a very narrow gap at Ep (of the order of 50 meV) for 
majority-spin and a much larger gap for minority spins. 
The spin moment, concentrated mainly at the V atoms, 
is an integer of exactly M — 2 fiB per unit formula, which 
is obvious by counting the occupied bands of the two spin 
directions. 

The exchange constants of VAs were calculated within 
GGA and adiabatic spin dynamics. The energy E(q) of 
frozen magnons is calculated as a function of the wavevec- 
tor q using the full-potential linearized augmented plane 
wave method (FLAPW) [3 which allows to evaluate 
the real-space exchange constants Jij by a Fourier trans- 
form of E{q). The procedure is similar to the one used 
by Halilov et al. |l9j. Using these exchange parameters 
in a Monte Carlo simulation of the corresponding classi- 
cal Heisenberg Hamiltonian E = — (1/2) JijMi ■ Mj 

(where Mj and Mj are the magnetic moments at sites i 
and j), we obtain a Curie temperature Tq = 820 K by the 
fourth-order cumulant crossing point. This resul t ag rees 
with the value of Tq — 830 K calculated in Ref. [20j us- 
ing a similar method. The high Curie point is well above 
room temperature making VAs a very promising candi- 
date for applications in spintronics. 

In order to take into account static correlations, we em- 
ploy the GGA+[/ method (using the values U = 2 eV, 
J = 0.9 eV typical for 3d transition metals [2lL H3|). 
Accordingly, we used the GGA equilibrium lattice pa- 
rameter, a = 5.69A, and a broadening 6 of about 15K, 
which allows the majority spin gap to be clearly resolved. 
For different lattice parameters (eg, the InAs parameter) 
our LDA results agrees with the ones previously reported 
Q|2^. The main difference between the GGA DOS (see 
Fig. [21 and GGA+J7 spectrum (not shown here) is, as 
expected, that within the GGA+U the occupied, local- 
ized majority e g states are shifted to even lower energy, 
while the unoccupied, minority e g states are shifted to 
higher energy. The semiconducting character does not 
change, since the e g and ti g bands remain separated for 
both spins; the majority-spin gap slightly increases but 
remains small. 

In order to investigate dynamic correlation effects in 
VAs we used a recently developed fully self-consistent in 
spin, charge and self-energy LSDA+DMFT scheme [23j |. 
This scheme uses for the LDA/GGA calculations the ex- 
act muffin-tin orbitals (EMTO) theory 0, and the full 
charge density technique. This method combines the ac- 
curacy of the full-potential method and the efficiency of 
the muffin-tin potential method [2J|. Correlation effects 
are treated in the framework of dynamical-mean-field 
theory with a spin-polarized T-matrix Fluctuation Ex- 
change type of quantum impurity solver 1251. |26|. The 
computational details are described in Refs. |23l 1271 128| . 

In our calculations, we considered the standard rep- 
resentation of the zinc-blende structure with an fee unit 
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FIG. 2: (color online) DOS of VAs within the GGA (dashed 
blue line) and GGA+DMFT (solid red line) for a temperature 
of T = 200A", U = 2 eV and J = 0.9 eV. Inset: Focus 
around Ef showing the semiconducting gap within the GGA. 
To ilustrate the minority-spin NQP states a ten times larger 
scale for the spin down channel is used. 



cell containing four atoms: V (0,0,0), As (1/4, 1/4, 1/4) 
and two vacant sites at (1/2, 1/2, 1/2) and (3/4, 3/4, 3/4) . 
The charge density is calculated by integrating the Green 
function along a complex contour up to the Fermi level. 
The multipole expansion of the charge density is cut off at 
Imax — 8, and the expansion of the Green function is cut 
off at l m ax = 3. The convergence of the calculations was 
checked up to a number of 1505 k vectors and several sets 
between 22 and 30 complex energy points. The PBE [3(J 
parameterization of the GGA exchange correlation po- 
tential was used. The one-particle, LSD A / GGA+DMFT 

Green function G a (k;E) is related to the LSDA/GGA 

-1 

and to 



Green function G°(k;E) 

the local (on-site) self-energy Y, a (E) via the Dyson equa- 
tion 



G-^k-E) =E + n-H a (k)-Il rT (E) 



(1) 



is the LSDA/GGA Hamiltonian, depen- 
011 the Bloch vector k, and the spin index 



H°(k) 
dent 

a 6 {Til}- /J is t ne chemical potential. The 
many-body effects beyond the LSDA/GGA are de- 
scribed by the multiorbital interacting Hamiltonian 

2 Ei{m,j) ^mm , ra"m" | C imJ C itn , (7 ,Ci m "V , Ci m « ff where 777, 

are local orbitals at site i; and c denote creation and 
destruction operators respectively. For the multiorbital 
Hamiltonian the on-site Coulomb interactions are 
expressed in terms of two parameters U and «/ 0- I n 
order to avoid "double counting" , the static part of the 
self-energy is subtracted, i. e. T, a (E) is replaced with 
T, a (E) — S CT (0), as S,t(0) is already included in the 
LSDA/GGA part of the Green's function. It has been 
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proven that this type of "metallic" double-counting is 
suitable for medium correlated d-electron systems [29(. 

The computational results for the GGA and 
GGA+DMFT densities of states are presented in Fig. 
The non-quasiparticle (NQP) states in the minority spin 
band are visible just above the Fermi level (inset), pre- 
dicted also by previous calculations jU • The weak 
spectral weight of NQP state is due to the fact that the 
Fermi level is close to the right edge of the minority spin 
gap, as discused for CrAs having a similar structure [28j • 
The origin of the NQP states is connected to the "spin- 
polaron" processes: the spin-down low-energy electron 
excitations, forbidden for the HMF in the one-particle 
picture, turn out to be possible as superpositions of spin- 
up electron excitations and virtual magnons || E| . The 
local spin moments at the V atoms do not change signif- 
icantly (less than 5%). 

However, in the case of VAs another correlation ef- 
fect appears: the small majority-spin gap at Ep closes, 
making the material half metallic. This prediction is the 
central result of our work. In order to investigate the 
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FIG. 3: (color online) Energy dependence of real and imagi- 
nary parts of the self-energy E CT (i?) for the tig orbitals. 

mechanism of the gap closure for the majority spin chan- 
nel, we look at the behavior of the electron self-energy. 
The quadratic form of the imaginary part of the majority- 
spin self-energy, ImEf (i?) ~ (E — Ep) 2 , visible in both 
Figs. |31 and indicates a Fermi liquid behavior, as op- 
posed to T,i(E) which shows a suppression around Ep 
due to the band gap, as well as a peculiar behavior for 
E > Ep related to the existence of NQP states. 0, 0. 
From the Dyson equation QJ, one can see that the real 
part of the self-energy, KeS a (E), causes a shift of the 
LDA energy levels. Therefore, due to the non-zero E^ 9 , 
the e g orbitals in the close vicinity of the Fermi level are 
pushed closer to Ep. This renormalization of the levels is 
connected to the large value of Re(<9E /8E)e f < 0, push- 
ing states in the vicinity of Ep more closer to Ep. This 
causes occupied levels to be renormalized to higher en- 
ergy and unoccupied levels to lower energy. Notice that 
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FIG. 4: (color online) Energy dependence of real and imagi- 
nary parts of the self-energy Yi a (E) for the e fl orbitals. 



this effect is completely opposite to the GGA+U results, 
discussed above. In addition to this shifting, the e g peak 
is broadened by correlations, its tail reaching over the 
Fermi level (Fig. [21 inset). Thus, our finite-temperature 
GGA+DMFT calculations demonstrate the closure of the 
narrow gap in the spin-up channel, which is produced by 
the correlation-induced Fermi- liquid renormalization and 
spectral broadening. At the same time, NQP states ap- 
pear for the minority spin channel just above Ep. 

The slope of the majority spin channel self-energy 
is almost a constant as a function of temperature: 
Re(dT,^/dE) BF ~ -0.4 between 200 K and 500 K. The 
physical content of this lies in the quasiparticle weight 
Z = (1 — <9ReE|/c9-E) _1 , measuring the overlap of the 
quasiparticle wave function with the original one-electron 
wave function, having the same quantum numbers. Our 
numerical results indicate that Z « 0.7 is quite tem- 
perature independent for small temperatures. Thus, we 
could expect this renormalisation to hold down to zero 
temperature. As a consequence the closure of the gap 
in the majority channel is a quantum effect, originating 
from the multiorbital nature of the local Coulomb inter- 
action (orbitals are squeezed towards Ep) rather than 
an effect of temperature. We have verified that a similar 
gap closure is obtained for larger values of U, namely, 
[7 = 4 and 6 eV, although the latter values should be 
taken with some caution in our FLEX calculation, which 
is in principle appropriate only in weak to intermediate 
coupling. As a general tendency, increasing U — J pro- 
duce a stronger Fermi-liquid renormalization in the ma- 
jority spin channel. We have also verified that the same 
effect is evidenced for J = Oel^. 

Recent mean-field (LDA+U) calculations by Anisimov 
et al. [3lJ yield a first-order, semiconductor to ferromag- 
netic metal transition as a function of doping in the 
FeSii-^Gea, alloy. In contrast, our calculation clearly 
shows that in VAs the closure of the semiconducting gap 
for majority spins cannot be captured by a static ap- 
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proach; dynamic correlation contributions in the multi- 
orbital model are required. 

In summary, we have investigated the electronic struc- 
ture and correlation effects in the zinc-blende alloy VAs. 
On the one hand, our density-functional theory calcula- 
tions within the GGA predicts this material to be a ferro- 
magnetic semiconductor with a tiny gap of about 50 meV 
in the majority-spin DOS. On the other hand, dynami- 
cal effects described by LDA+DMFT destroy the narrow 
band gap and turn the material into an half-metallic fer- 
romagnet. According to our results the closure of the 
band gap is due to the multiorbital nature of the local 
Coulomb interaction, and can be described as a strong 
correlation-induced Fermi liquid-like renormalization of 
majority-spin states accompanied by a lifetime broad- 
ening. At the same time, in the minority-spin channel 
non-quasiparticle states appear just above Ep. We stress 
that these results are a consequence of the interaction 
of spin-up electrons with spin- flip excitations. Neither 
the closure of the gap nor the non-quasiparticle states 
can be obtained in the static LDA+J7 approximation. 
Our LDA/GGA calculation supplemented by a Monte 
Carlo simulation also predicts a high Curie temperature 
of 820 K, which makes this material of interest for tech- 
nological applications. We expect T c to be not much 
affected by dynamical correlation, for the same reason 
for which the effective exchange interaction parameters 
are not affected, as it was demonstrated in recent works 

The revealed half-metallic (instead of semiconduct- 
ing) behavior has important consequences in the poten- 
tial applications of VAs in spintronics. In contrast to 
all-semiconductor-based spin-injection devices Q which 
avoid the resistivity mismatch problem, half-metals can 
be applied in giant magnetoresistance or, if interface 
states are eliminated |23, m tunneling magnetoresis- 
tance. As our calculations shows, in the prediction of 
new spintronic materials, correlation effects play a deci- 
sive role. While in some materials these are detrimen- 
tal for half-metallicity due to the introduction of spec- 
trum in the minority spin gap [22], the present case is 
an example in which correlations turn out to be favor- 
able for a high spin polarization. The metallic nature of 
the majority spin channel would be visible in resistivity 
measurements. Therefore, the experimental realization 
of zinc-blende VAs would provide a test of our predic- 
tion. Further research should address the issue of the 
stability of the half-metallic ferromagnetic state in a zinc 
blende structure. Some work in this direction has been 
already carried out [3J, |35j . DMFT total-energy calcula- 
tions along the lines described in Ref. [2a| are in progress. 
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